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The expansion of the strongly interacting quark gluon plasma (sQGP) created in Au+Au colli- 
sions at RHIC can be described by hydrodynamical models. Hadrons are created after a freeze- 
out, thus their distribution describes the final state of the evolution. The earlier stages can be 
analyzed via penetrating probes like photon observables. These were measured in 2010 and 201 1 
by the PHENIX experiment. Here we analyze an analytic, 1+3 dimensional perfect relativistic 
hydrodynamic solution and calculate hadron and photon observables, such as transverse momen- 
tum spectra, elliptic flow and correlation (HBT) radii. We find that our model is not incompatible 
with the data, not even with the direct photon elliptic flow. From fitting the data, we find that early 
temperatures of the sQGP were well above the quark-hadron transition temperature, in the hottest 
point, the center of the fireball the temperature may have reached 507±12 MeV. The equation 
of state of this quark matter can be described by an average sound speed of 0.36±0.02. We also 
predict a photon source that is significantly larger in the out direction than in the side direction. 
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1. A hydrodynamic model 

It is a well established picture that a strongly interacting quark gluon plasma is created 
in relativistic Au+Au collisions of the Relativistic Heavy Ion Collider, and the evolution of this 
medium can be described by perfect hydrodynamics. The equations of hydrodynamics can be 
solved numerically, which has the advantage of having arbitrary initial conditions. It is also possible 
to find analytic solutions with realistic properties. However, very few truly 1+3 dimensional (and 
not spherically symmetric), relativistic models were compared to data yet. In this paper we analyze 
the relativistic, ellipsoidally symmetric model of Ref. [^. Hadronic observables were calculated in 
Ref. [^, while photonic observables in Ref. [||]. 

The picture used in hydro models is that the pre freeze-out (FO) medium is described by hydro- 
dynamics, and the post freeze-out medium is that of observed hadrons. The hadronic observables 
can be extracted from the solution via the phase-space distribution at the FO. This will correspond 
to the hadronic final state or source distribution S{x,p). See details about this topic in Ref. |^]. 
It is important to see that the same final state can be achieved with different equations of state or 
initial conditions However, as discussed below, the source function of photons is sensitive to 
the whole time evolution, thus both to initial conditions and equation of state as well. 

For the direct photon calculations [^|, the key assumption is, that even though direct photons 
may not be thermalized in the strongly interacting plasma (as their mean free path may be on the 
order of the size of the fireball), but the radiation itself is thermal. Thus the phase-space distribution 
of the photons is characterized by the temperature of the medium (at a given space-time cell), while 
the expansion of the fireball also effects the observed spectrum. This is a macroscopic model, and 
in the following we will calculate photon observables from it and compare to RHIC data. The most 
important assumption is, that the spectrum of direct photons is thermal because macroscopically, 
the photon radiation is thermal. 

The analyzed solution [Q] assumes self-similarity and ellipsoidal symmetry, as described also 
in Refs. Q]. The ellipsoidal symmetry means that at a given proper time the thermodynamical 
quantities are constant on the surface of expanding ellipsoids. The ellipsoids are given by constant 
values of the scale variable s, defined as 

rl r^y r] 

■"r^-r^, (1.1) 



W y2f2 22^2^ 

where the constants X, F, and Z describe the expansion rate of the fireball in the three spatial 
directions. Spatial coordinates are r,^, Yy, and r^. The velocity-field is described by an izotropic 
Hubble-type expansion: 

ui'{x) = — (1.2) 



where x means the four-vector (t,rx,ry,r^) and T = Vx^ is the proper-time coordinate. 
The temperature distribution T{x) is given as 

T{x) = To(^y^\xp(^], (1.3) 
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where T is the proper time, s is the above scahng variable, while Tq = r| v=o,T=To> and To is an 
arbitrary proper time, but practically we choose it to be the time of the freeze-out, thus Tq is the 
central freeze-out temperature. Parameter b is proportional to the temperature gradient, i.e. if the 
fireball is the hottest in the center, then b <0. If there is a conserved charge in the system e.g. 
the baryon number density, then charge number density n{x) can be utilized in the solution. As 
described in Refs. ^], such a number density can be introduced as 

n(,)=nJ^)\xp(-^-^). (1.4) 



For the momentum distribution of direct photons, this will not be needed, as the only the tem- 
perature of the medium (the strongly interacting plasma) governs the creation of photons, not the 
density (which however plays an important role also in the case of hadron creation). The equation 
of state (EoS) we use here is e = Kp, with e being the energy density and p the pressure. Here 
K = (one over speed of sound squared) is the main parameter describing the EoS. 

From the above hydrodynamic quantities, source functions can be created. For bosonic hadrons, 
it takes the following form : 

W4 p^d'Z^{x)H{T)dr 

S{x,p)d x = J^— — -, (1.5) 

n{x)exp[p^uf^{x)/T[x)) — 1 

where ^ = g/{2n)^ (with g being the degeneracy factor), H{t) is the proper-time probability 
distribution of the FO. It is assumed to be a 5 function or a narrow Gaussian centered at the freeze- 
out proper-time To. Furthermore, n{x)/T{x) = InM(x) is the fugacity factor and d^l.^{x)p^ is the 
Cooper-Frye factor [^] describing the flux of the particles, and d^'L^{x) is the vector-measure of 
the FO hyper-surface. Here the source distribution is normalized such as / S{x,p)d^xd^p/E = N, 
i.e. one gets the total number of particles N (using c=\,%=\ units). 
For the source function of photon creation we have [H] : 



p^d'^l.^{x)dt _ p^ 



s{x,p)d\ = — " r'r,, — = — -, — , — d'^x (1.6) 

exp [p^u'^{x)/T{x)) - 1 exp [p^u^{x) /T [x)) - 1 

where p^d^L^ is again the Cooper-Frye factor of the emmission hypersurfaces. Similarly to 
Ref. [||] we assume that the hyper-surfaces are parallel to u^, thus d^'L^{x) = u^d^x. This yields 
then p^u^ which is the energy of the photon in the co-moving system. The photon creation is the 
assumed to happen from an initial time ti until a point sufficiently near the freeze-out. 

Experimental observables can then be calculated from the source function, using a second 
order saddlepoint approximation. In this approximation the point of maximal emissivity is 

Pi 

ro.i = Pit — Jot i = x,y,z (1.7) 
while the widths of the particle emitting source are 

/?? = PrToj(-] ,for/ = x,y,z (1.8) 

where we introduced the auxiliary quantities 

Pi = (1-9) 

K-'i-Kb/R^ 

where again k: = cj^ is describing the EoS, and /?, = X,F,Z for / = x,y,z, respectively. 
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2. Calculated observables 

The invariant one-particle momentum distribution is defined as Ni{p) = J S{x,p)d'^x. It de- 
pends on tiie three-momentum p = {px,Py,Pz)- We will introduce the {pt,(p,p^) cylindrical coor- 
dinates (z being the beam direction) and use the longitudinal rapidity y (for which Edy = dp^ is 
true). As usual, we will restrict our calculations to j = (note that in this case E = pt is true for 
photons). Our calculated quantities will then be the elliptic flow V2, and the transverse momentum 
distribution N\ (pt). These can be calculated from Ni (p) as 

Ni{pr)= A^i(;^)|,.o^'P (2-1) 

''^''^ = md^) ^'-'^ 

We also calculated Bose-Einstein correlation radii from our model. As usual, the two-particle 
correlation function for identical particles can be calculated from the single particle source function 

S{x,p) as 



C2{q) = l+X 



S{q) 



S{0) 



(2.3) 



where q is the momentum difference of the two particles and S{q) is the Fourier-transformed of 
the source S{x,p) in the variable x, and the dependence on the momentum p is not noted in the 
formulas. This correlation function has, as usual, a shape with a peak, the width of which is 
characterized by the HBT radii /?out, ^side and /?out- We calculated these radii for different average 
momenta p. 

Here we do not detail the analytic result of these calculation, but will show a comparision of the 
model to the data in the next section. The detailed results are given in Refs. [^, Q]. It is important to 
note however, that in the final formulas, we use transverse expansion (m,) and eccentricit (e) instead 
of X and y direction expansion rates X and Y: 

-^ = -{^ + ^],e = -. ^. (2.4) 



3. Comparison to the measured direct photon spectrum 

The freeze-out parameters were determined from hadronic fits in Ref. These properties 
include the expansion rates, the freeze-out proper-time and freeze-out temperature (in the center of 
the fireball), as shown in Table |l|. When describing direct photon data, we used the parameters of 
the hadronic fit and left only the remaining as free parameters [Q]. The free parameters will be K 
(the equation of state parameter) and ti, the initial time of the evolution. 

We compared the above results to PHENIX hadron and photon data of 200 GeV Au-i-Au col- 
lisions. We fitted our above formulas to PHENIX invariant transverse momentum distributions of 
Ref. [^, HBT radii of Ref. ^ and elliptic flow data of Ref. [Q]. We used direct photon data also 
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AJ, and HBT 


pllintir flow 








0-30% cent. 


0-30% cent. 


0-92% cent 












r^ptitral PO tpmnpratiirp 


7n TMeVl 


199±3 


204 ±7 


204 MeV Cn xed') 


Eccentricity 


£ 


0.80±0.02 


0.34±0.03 


0.34 (fixed) 


Transverse expansion 


uj/b 


-0.84±0.08 


-0.34±0.01 


-0.34 (fixed) 


FO proper-time 


To [fm/c] 


7.7±0.1 




7.7 (fixed) 


Longitudinal expansion 




-1.6±0.3 




-1.6 (fixed) 


Equation of State 


K 






7.9 ±0.7 


Initial time 


tj [fm/c] 






0-0.7 fm/c 


Fit quality 


Degrees of freedom 


NDF 


41 


34 


3 


Chisquare 




24 


66 


7 


Confidence level 




98% 


0.1% 


7.2% 



Table 1: Parameters of the solution, describing the expanding sQGP. The first five were determined from 
hadronic fits the remaining from direct photon data [Q]. See details in the text and in these references. 
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Figure 1: Fits to invariant momentum distribution of pions (top left), HBT radii |Q| (top right) and 
elliptic flow (bottom). See the obtained parameters in Table In the middle plot the lower curve is the 
"aiXaRout and R^e- 
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Figure 2: Fit to direct photon invariant transverse momentum data [|10p (top left), comparision to elliptic 
flow data [|llj (top right) and direct photon HBT predictions (bottom). See the model parameters in Table [l| 

from PHENIX |ll|]. Results are shown in Figs. |l| and ^, while the model parameters are detailed 
in Table. |l|. 

Let us discuss the results of the photon fits! The equation of state result is K" = 7.9 ±0.7, or 
alternatively, using K = 



0.36 ±0.02,,„,± 0.04 



syst 



(3.1) 



which is in nice agreement with both lattice QCD calculations [ |12| ] and experimental results from 



hadronic data [ ]13| , |14| ]. This represents an average EoS as it may vary with temperature. There may 
be solutions with a k{T) function [p^], but for the sake of simplicity we assumed here an average, 
fixed K. As detailed in Ref. [Q], we determined an "interval of acceptability" for f,. The maximum 
value for ti within 95% probability is 0.7 fm/c. This can then be used to determine a lower bound 



for the initial temperature, using the eq. |1.3[ Thus the initial temperature of the fireball (in its 
center) is: 

Ti > 507 ± n,rat ± 90,,.,,MeV (3.2) 
at 0.7 fm/c. This is in accordance with other hydro models as those values are in the 300 — 600 MeV 



interval [10]. Note that a systematic uncertainty was determined by using a prefactor of {T /Tq)^, 
with A'^ = 0, 1,2,3. This factor arises if the photon creation can be described by a microscopic 
process, as detailed in Ref. This causes only a minor change in the resulting spectrum, as it 
is dominated by the exponential factors in it. However, the equation of state parameter K changes 
from 7.9 to 6.5 as we increase the exponent in the prefactor. This gives a systematic uncertainty to 
our parameters. 
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A measurement of direct photon elliptic flow was also performed recently at PHENIX [pl|]. 
Using the previously determined fit parameters we can calculate the elliptic flow of direct photons 
in Au+Au collisions at RHIC. Due to the low number of points in the desired range, a fit could not 
be performed here, but we used the average value e in case of the two fits of Ref. [^]. The resulting 
curve, where the value e = 0.59 was used, is shown on Fig. ^ 

In case of hadronic HBT, correlation radii in the side and out directions are almost equal, as for 
the hadronic transition is of cross-over type (i.e. the transition time is short), see details in Ref. [^]. 
However, in case of photons, the creation spans the whole evolution of the fireball, thus /?out will 
be significantly larger than /?side- Indeed this was observed in our model, as shown on Fig. ^ 

To summarize, we find that thermal radiation is consistent with direct photon data, and our 
result on the equation of state is = 0.36 it 0.02^^^^ it 0.04^3,^^. We set a lower bound on the initial 
temperature of the sQGP to 507 ± ll^tat ± 90^^,^? MeV at 0.7 fm/c. We also find that the thermal 
photon elliptic flow from this mode is not incompatible with measurements. We also predicted 
photon HBT radii from the model, and discovered a significantly larger /?out than /?side- 
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